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Measurement of Neutron Background at the Pyhasalmi mine for CUPP 
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J. Peltoniemi b , T. Keranen b 

a Institute for Nuclear Research, Russian Academy of Sciences, Moscow, 117312, Russia 
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A natural neutron flux is one of significant kind of background in high-sensitive underground experiments. 
Therefore, when scheduling a delicate underground measurements one needs to measure neutron background. 
Deep underground the most significant source of neutrons are the U-Th natural radioactive chains giving a fission 
spectrum with the temperature of 2-3 MeV. Another source is the U-Th (a,n) reactions on light nuclei of mine 
rock giving neutrons with different spectra in the 1-15 MeV energy region. Normal basalt mine rocks contain 
10~ 6 g/g of 238 U and less. Deep underground those rocks produce natural neutron fluxes of 10 -7 -10 -6 cm _2 s _1 
above 1 MeV. To measure such a background one needs a special techniques. 

In the Institute for Nuclear Research, Moscow, the neutron spectrometer was developed and built which is 
sensitive to such a low neutron fluxes. At the end of 2001 the collection of neutron data at the Pyhasalmi 
mine was started for the CUPP project. During 2002 the background and rough energy spectra of neutron at 
underground levels 410, 660, 990 and 1410 m were measured. The result of the measurement of the neutron 
background at different levels of the Pyhasalmi mine is presented and discussed. Data analysis is performed in 
different energy ranges from thermal neutrons up to 25 MeV and above. 



1. Introduction 

When developing underground physics labora- 
tories it is need to take into account the fac- 
tors which contribute in background environ- 
ment. One such factor is neutrons which are pro- 
duced by natural radioactive U-Th chains in sur- 
rounding rocks. At depths which are high enough 
(2000 m.w.e. and more) the intensity of cosmic 
ray muons decreases down to ~10~ 7 cm _2 s -1 
which is comparable to the background neutron 
flux from crystalline rocks. Thus, at high depths 
a relative contribution of neutrons to background 
of physics experiments may appear to be con- 
siderable. Estimation of this contribution into 
background can be based on calculation of flux 
using the measured concentrations of U-Th and 
chemical composition of surrounding mine rocks. 
Direct measurement of the spectrum and flux of 
background neutrons is important too. However, 
a low neutron flux is very difficult to be measured. 

'Corresponding author: jna@al20.inr.troitsk.ru 



Several research groups have investigated the 
neutron background at different underground lab- 
oratories |1I2I3| . Some of them used 6 Li-dopped 
liquid scintillator technique |2] and others used 
in addition a pulse shape discrimination tech- 
nique 1 . The fluxes of ~10 -6 cm~ 2 s~ 1 in 
the energy range >2 MeV were measured there. 
Sometimes a radiochemical technique was used 
there |3j. 

In the Institute for Nuclear Research, Moscow, 
the neutron spectrometer was developed and built 
which is sensitive to such a low neutron fluxes. 
The measurements of fast neutron flux in several 
underground rooms of Baksan Neutrino Observa- 
tory (BNO), SAGE and DULB jSJ, have been 
performed with the spectrometer. These rooms 
are located under Mt. Andyrchy (Northen Cau- 
casus Mountains, Russia) in a tunnel that pen- 
etrates 4.5 km into the mountain, at a depth of 
4700 and 4900 meters of water equivalent accord- 
ingly. The fluxes of ^10~ 7 cm _2 s _1 in the en- 
ergy range >1.5 MeV were measured there suc- 
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Figure 1. General view of the detector 



cessfully. 

The series of neutron flux measurement were 
carried out according to agreement for a joint 
Russian-Finnish scientific research programme 
in 2000-2002 between the Institute for Nuclear 
Research of the Russian Academy of Sciences 
(project SAGE) and The Oulu Southern Insti- 
tute of the University of Oulu (project CUPP, 
http://www.cupp.oulu.fi) at underground facili- 
ties of the Pyhasalmi mine, town Pyhajarvi, Fin- 
land. We present the results here. 



2. Detector 

2.1. Design 

The detection part of the spectrometer (detec- 
tor) consists of a liquid organic scintillator viewed 
by three photomultipliers (PMT) and nineteen 
proportional counters filled with 3 He (neutron 
counters — NC) and distributed uniformly over 
the scintillator. Figure Q shows a general view of 
the detector. It was described initially in [7|. 

The housing of the detector is made of stain- 
less steel. The space of cylinder tank (diameter 
of 360x360 mm long) available for the scintilla- 



tor is about 30 1. To enhance light collection the 
wall and bottom of the scintillator tank are coated 
with a fluoroplastic. Each counter in the tank is 
coated with aluminium foil too. The illuminator 
is made of organic glass and is covered by fluo- 
roplastic plate with windows intended to improve 
the light collection and to secure the phototubes. 
The external surface of the tank is 6267 cm 2 . 

2.2. Principle of operation 

Fast neutrons (>1 MeV) entering the scintilla- 
tor are decelerated down to a thermal energy and 
diffuse in the detector until they are either cap- 
tured in a neutron counter or captured by scin- 
tillator protons or leave the detector. The inten- 
sity of light flash combined from recoil protons 
which are produced during neutron thermaliza- 
tion is on average proportional to the initial neu- 
tron energy. A portion of thermalized neutrons 
are captured by 3 He nuclei in the neutron coun- 
ters, which emit charged particles via the reaction 
3 He(n,p)t, E p =574 keV, E t =191 keV. Thus, a sig- 
nature of such event is a light flash in the scintilla- 
tor followed by capture in the neutron counter af- 
ter some time of delay. This delay is conditioned 
by the mean lifetime of thermalized neutrons in- 
side the detector and is determined mainly by the 
detector design. For this particular detector T 1 /^ 
is ^50 /is. This technique allows one to suppress 
the natural background of 7-rays significantly, by 
several orders of magnitude. 

2.3. Acqusition system 

The functional diagram of the data acquisition 
system is shown on Figure In order to simplify 
the apparatus design the signals from all PMTs 
and NCs are mixed into two independent channels 
called "PMT channel" and "NC channel" respec- 
tively. The circuit design of the low-level signal 
electronics are selected taking into account an op- 
timal signal-to-noise ratio. To reduce the noise, 
the main amplification of initial signals is accom- 
plished with the preamplifiers before joining. 

At the PMT channel a negative signal from 
anode of the PMTs is entered to the input of 
preamplifier (PA). A continuously adjusted PA 
gain ensures the operation of all three PMTs from 
a single high-voltage power supply. The PMT 
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Figure 2. The functional diagram of a data ac- 
quisition system 



signal from PA output is multiplexed by a fast 
inverting linear mixer and then is distributed on 
two directions. One branch entered directly to 
the input of first channel of digital oscilloscope 
(DO). The other one is entered to the input of 
combined unit of amplifier and single-channel an- 
alyzer (SCA). The positive TTL-specified signal 
from SCA output triggers the DO through an ex- 
ternal trigger input at the " Spectrum Acquisition 
Mode" (see below) which is used in calibration 
procedures, control of the system operation sta- 
bility and background rate monitoring during ac- 
quisition. 

In the NC channel, a high voltage of positive 
polarity from a single source is applied through 
high-voltage isolating resistors to the anodes of 
helium counters. The anode generates signal 
of negative polarity which is fed through high- 
voltage separating capacitors to the input of PA. 
The spread of the counter gas amplification fac- 
tor is compensated by adjusting the PA gain dur- 
ing calibration. For convenient channel tuning 
each of counter PA (similarly to the PMT chan- 
nel) can be switched off independently by turning 
off the PA supply voltages. The signal of neu- 
tron counter PA outputs are multiplexed in a sin- 
gle output fast inverting linear mixer and subse- 



quently enter directly to the input of the second 
channel of DO. A signal from NC channel trig- 
gers the data acquisition system at the "Pulse 
Acquisition Mode" which is usually used for real 
background measurements. 

The full waveform of event in the PMT and 
NC channels is recorded independently by means 
of the two-channel PC/AT-interfaced digital os- 
cilloscope LA-nlOM5 manufactured by "Rudnev- 
Shilyaev" Co., Moscow, Russia. It's performance 
data are shown in Tabled The frame is recorded 
in selected time interval including periods before 
and after the trigger which are called prehistory 
and history accordingly. This time interval can 
be adjusted in the acquisition algorithm over the 
wide range (2 (j,s to 2 ms). Usually the trigger, 
that is a capture in 3 He counter, is put in the 
middle of scope frame (—160 to 160 /is). Thus, 
the prompt PMT event being correlated with the 
3 He count should appear before it in the interval 

— 160 to /xs and uncorrelated one may appear 
elsewhere with the flat probability over the whole 
frame. 

The acquisition system operates under control 
of code NOSC especially designed for this pur- 
pose. It is written in C and is executed in real- 
time mode MS-DOS 6.22. Mainly, it is a virtual 
oscilloscope recording frame by frame in "Pulse 
Acquisition Mode" . In addition, it is able to col- 
lect a pulse height histogram in "Spectrum Ac- 
quisition Mode" . 

2.4. Sensitive elements 

The scintillator "Protva" is developed and 
manufactured by M. Konoplya, IHEP, Protvino, 
Russia. It is similar to NE-213 and has the den- 
sity of 0.84 g/cm 3 , the hydrogen concentration 
of 136 g/kg, the flash duration of 3 ns, the light 
yield at least 40% of antracene and the ignition 
temperature of 80° C. The scintillator is viewed 
with three phototubes PMT Model 173, manu- 
factured by "Ekran", Novosibirsk, Russia. The 
photocathodc of PMT is 160 mm in diameter, so 
three tubes cover about 70% of top side of the 
tank. The PMT Model 173 isn't designed for fast 
timing applications, so the rise time is rather poor 

— about 100 ns. The preamplifiers of signals from 
PMTs and the signal mixer are placed at the top 
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Table 1 



Main performance data of the digital oscilloscope LA-nlOM5 
Parameters are presented as for 2-channel (1-channel) mode. 



Independent A/D converter 


2 


Resolution of A/D conversion 


8 bit 


Sample frequency range 


3.125 kHz-50 MHz (6.104 kHz-100 MHz) 


Internal RAM 


256 kB, 128 kB per channel 


Input sensitivity, V 


±1, ±0.5, ±0.2, ±0.1 


Input impedance 


1 MQ and 15 pF 


S/N ratio 


47 dB 



of the detector and are covered with an optically 
isolated cap. 

Neutron counters are commercially available 
SNM-18, Russia. Nineteen counters (diameter of 
30x300 mm long) filled with a 3 He+4%Ar mix- 
ture at a pressure of 4 atm are mounted in sealed 
wells in the internal space of the tank. The width 
of current signal is about l-=-4 fia. The anode 
contacts of the counters are terminated under the 
bottom of the tank. The preamplifiers of NCs and 
the signal mixer are placed at the bottom of the 
detector and are also covered with a cap. 

2.5. Efficiency 

The efficiency was simulated in two stages. 
First, the history of the neutron with the en- 
ergy Eq initially incident on the surface of the 
detector was traced individually taking into ac- 
count all kinds of an energy loss. A special code 
DGEOM was developed there to trace the neu- 
tron. It traces the neutron down to 50 keV fixing 
those recoil protons which are able to produce a 
detectable light flash. It fixes those events where 
the neutron leaves the detector out also. As a re- 
sult one finds a geometrical efficiency g(Eo), i.e. 
the probability of the neutron to be trapped in 
the detector depending on initial energy Eq ■ The 
DGEOM does a correction of probability of inci- 
dent neutrons by a cosine of falling angle, so one 
can use a geometrical surface of the detector in a 
data analysis. 

Second, a standard code MORSE was used to 
simulate the efficiency Sth of capture of neutron 
in the detector. It's assumed that neutrons are 
burned in the detector uniformly with the energy 
50 keV — a final energy of previous stage. It is 
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Figure 3. Efficiency plot 



resulted in e t h = (20 ± 1)%. It was checked with 
Pu-Be source and found to be agreed within 10% 
uncertainty. 

Last, the efficiency of the detector to trap the 
neutron and to measure it's energy is found as 

e{E ) = e th x (Eq). (1) 

It is shown on the Figure 01 

3. Calibration of energy scale and time of 
delay 

3.1. NC channel 

The Pu-Be neutron source with intensity of 
2 x 10 3 s _1 was used to irradiate the detector 
keeping in mind an unification of individual spec- 
trum of each 3 He proportional counter. The pulse 
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Figure 4. The pulse height spectrum of mixed 
helium counters events when irradiated with Pu- 
Bc neutron source 



height spectrum of mixed nineteen counters is 
shown on Figure 01 Two peaks are clearly ex- 
pressed there: first one is located near 191 keV 
of tritium energy and second one is located near 
760 keV of both tritium and proton energies. The 
neutron window of acceptance is choosen to be 
170-870 keV, that are 35-180 channels of ADC in 
±100 mV range. The window accepts 90% of to- 
tal captured neutrons. The 760 keV peak is found 
at 157 ± 2 channel of ADC, so the weight of one 
ADC channel in ±100 mV range is 4.84 keV at 
1270 V of high voltage supplied. The uncertainty 
of the scale is estimated to be ±1%. 

3.2. PMT channel 

A 60 Co 7-source (1.17 and 1.33 MeV lines) was 
used to calibrate the PMT channel. Due to low Z 
the 7-ray absorption ability of the scintillator is 
rather weak (di/ 2 ~ 12 cm -1 for E g = 1 MeV). 
In the frame of particalur design of the detector 
one can't observe the peak of full absorption of 
60 Co. Therefore, a Compton edge location was 
used to calibrate the PMT channel. The detecor 
response on 7-rays was simulated there, and one 
found that the peak of Compton edge should ap- 
pear at 1.06 MeV. Figure El shows a pulse height 



3500 




Figure 5. The pulse height spectrum of mixed 
PMT events when irradiated with 60 Co 7-source 



spectrum of PMT signals when irradiated with 
60 Co. It resulted in the weight of one ADC chan- 
nel to be equal to 26.5 keV at 950 V of high volt- 
age supplied. One should note here that there is 
some nonuniformity of light collection. The same 
light flashed both at the top and at the bottom of 
the tank will result in different PMT signals. The 
difference was measured to be 20%, so it is the 
main contribution in the uncertainty. Thus, the 
uncertainty of an electron equivalent energy scale 
of the PMT channel is estimated to be ±10%. 

The relation between electron and neutron en- 
ergy scales was obtained as a secondary result of 
the efficiency simulation (see Figure [SJl keeping in 
mind a specific light yield of proton in the scin- 
tillator. This light yield is provided in [8] and is 
expressed as 

E 1 = 0.95E P - 8(1 - exp(-0.1£;0- 9 )), (2) 

where E 1 and E p are in MeV. It results in nega- 
tive light yield for E p < 0.2 MeV, so we use sim- 
ple quadratic form for the E p < 0.5 MeV energy 
range. It was found that the light flash produced 
by recoil protons scattered by neutron with E n is 
on average equivalent to the light flash emitted 
by electron with E e , which is expressed as 

E e = 0.134£^- 421 , (3) 
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Figure 6. The electron v. neutron energy scale 
for the NE-213 scintillator 



where E e and E n are in MeV. 

Finally, the energy of neutrons expressed in 
ADC channels and electron equivalent MeV is 
presented in Tabled 

3.3. Correlated signals 

The correlated events from the detector irradi- 
ated with Pu-Be neutron source are used to find 
a specific lifetime of thermalized neutrons before 
capture in helium counters. The delay time distri- 
bution for the neutron source is presented on Fig- 
ure [7| The best fit of lifetime is 2\ / 2 = 50 ± 2 //s. 

4. Data acquisition 

We started the collection of neutron data at 
the Pyhasalmi mine at the end of 2001. Dur- 
ing 2002 we were measuring the background and 
rough spectra of neutron at underground levels 
410, 660, 990 and 1410 m. In the TableElone can 
find history of measurements. The files have be- 
ing recorded 55 min in "Pulse Mode" and 5 min 
in "Spectrum Mode" each hour. 

There are two modes of measurement — LOW 
and HIGH (see Table |2J which differs mainly by 
range of PMT ADC (±100 or ±200 mV). The 
LOW mode is used to find the neutron flux in 
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Figure 7. The time of delay distribution of fast 
neutron events (Pu-Be source) 



the 0-12 MeV range as well as the HIGH one is 
used for >12 MeV energies. 

The passive shield made from standard lead 
sheets (1000 x 500 x 10 mm) in four layers was 
used there. It is used to suppress by several times 
a background 7-rays producing accidental coin- 
cidences with a negligible distorting of the neu- 
tron spectrum. The 7-ray counting rate under the 
shield varied in 400-600 s _1 range above 100 keV 
over all levels of the mine. 

5. Events characterization 

5.1. Normal events 

Normal events are reasoned by particles pro- 
ducing both a light flash in the scintillator and 
current signal in the proportional counter. One 
significant criterium is a shape of pulses produced 
by the preamplifiers. Normal PMT event should 
have the rise time to be ~0.1 fis which is reasoned 
by the time of electron collection in the photo- 
multiplier. The fall time is reasoned only due to 
discharge of correspondent input and output RC 
circuits of the preamplifier and has width ^30 fis. 
Normal NC event has, in contrast, the rise time 
l-j-4 /is due to internal properties of the propor- 
tional counter. Again, the fall time is reasoned 
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Table 2 

Intervals of neutron energy expressed in different units for two modes of measurement 

Mode Units Neutron energy interval 

MeV, neutron scale 0-1.5 1.6-3 3-6 6-12 12-25 >25 

LOW MeV, electron scale 0-0.25 0.26-0.62 0.63-1.65 1.66-4.43 1 1 

ADC chan., ±100 mV 0-9 10-24 25-63 64-168 

HIGH MeV, electron scale 4.4-12.6 >12.6 

ADC chan., ±200 mV 84-255 256 



Table 3 uo 



Chronology of the measurements 



"LOW" 


energy 


range means 


1.5-12 MeV, 


"HIGH" 


means 


>12 MeV. 




Level 


Mode 


Data files 


Date 


1410 


LOW 


J112212 


21 Dec-03 Jan 






J201032 


03 Jan-11 Jan 


440 


LOW 


J202130 


13 Feb-26 Feb 






J202270 


27 Feb-06 Mar 




HIGH 


J203060 


06 Mar-25 Mar 






J203270 


27 Mar-15 Apr 






J204150 


15 Apr-06 May 


990 


LOW 


J205070 


07 May-27 May 




HIGH 


J206110 


11 Jun-02 Jul 


660 


LOW 


J207020 


02 Jul-29 Jul 




HIGH 


J207290 


29 Jul-13 Aug 


1410 


HIGH 


J209131 


13 Sep-27 Sep 



only due to discharge of the RC circuits at the 
input and output of preamplifier and has width 
^60 /xs. 

Another significant criterium is a relative ap- 
pearance of signals in the frame of the scope. One 
example of normal event is shown on Figure |HJ It 
is a good candidate to be a neutron thermalized 
in the scintillator with prompt PMT signal and 
captured in one of the 3 He counters. Such event 
is most probably a correlated coincidence. An- 
other source of such events is a decay of Bi and 
Po radioactive isotopes, which can take place in 
a helium counter wall. It has been considered as 
main possible source of an internal background. 
It imitates an actual neutron event when beta 
decay of 214 Bi firing the scintillator is followed 
by delayed (Tj/2 = 164 /is) a-decay of 214 Po in 
helium counters. The delay time distribution of 
correlated events was measured in special water 
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Figure 8. Normal event with the PMT signal be- 
fore. Upper line — PMT, lower — NC 



shield suppresing the neutron background com- 
pletely [5]. It was found that the T X j 2 is very 
close to 164 /xs confirming thus the origin of the 
detector internal background. 

There is some probability for the event to be 
completely uncorrelated, or an accidental coin- 
cidence. The source of accidental coincidence 
is an independent light flash, produced by a 
background 7-ray with simultaneous capture of 
a background thermal neutron or a-particle in 
helium counter. For those events the PMT sig- 
nal may appear before as well as after the helium 
counter signal with the same probability. An ex- 
ample of such event is shown on Figure The 
helium counter signal is reasoned most probably 
by uncorrelated 214 Po decay here because of the 
amplitude doesn't fit the neutron window of ac- 
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Figure 9. Normal event with the PMT signal af- 
ter. Upper line — PMT, lower — NC 



ceptance (that is overflowed). A passive shield 
suppresses those events significantly. 

Normal event may appear also when a back- 
ground 7-ray produces an electron that fires both 
the scintillator and the helium counter at the 
same time. It results in the frame recorded where 
the PMT signal is located at the same time po- 
sition as the NC signal is. An example of such 
event is shown on Figure ITU1 

5.2. Breakdown, Discharge and Noise 
events 

Sometime a breakdown event in helium counter 
may appear there. An example of such event is 
shown on Figure^] One can easily note that the 
rise time of the event is rather fast in comparison 
to normal event. This is due to that a breakdown 
is a discharge appearing outside the counter just, 
for example, through water traces on the surface 
of counter. Due to variation of humidity in the 
mine air sometimes one observes high breakdown 
rate for a long time (several hours or more). 

An origin of discharge events in helium counter 
is most probably a discharge through air outside 
the counter. Due to low drift velocity it is rather 
slow one (see Figure l*H?ll . 

Noise events may appear due to high frequency 









■ 

\ 











, — I I I [ 1 

-150 -100 -50 50 100 150 



Delay, llsec 



Figure 10. Normal simultaneous event. Upper 
line — PMT, lower — NC 



current induced by power turning off/on or some 
other reasons. An example of the noise event is 
shown on Figure ITS1 

5.3. Discrimination 

It's easy to discriminate completely breakdown, 
discharge and noise events due to a very specific 
pulse shape of the signals. It's easy also to reduce 
a pure accidental coincidence when the PMT sig- 
nal appears after the NC one. This kind of dis- 
crimination is applied to data. 

In contrast, for the particular design of the de- 
tector it is practically impossible to distinguish 
a different sources of the individual correlated 
events. Nevertheless, it becomes quite possible by 
means of statistical method when one collect dif- 
ferent events at the same time of single measure- 
ment. So, a pure fast neutron signal is obtained 
as a remainder of the correlated data substracted 
with the accidental and Bi-Po background. 

6. Data analyzis 

An immediate result of measurement is a file 
containing the frames recorded. The analysis of 
data is proceeded in four stages. 

First, one analyses periods of the event rate 
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Figure 11. NC breakdown event. Upper line 
PMT, lower — NC 



Figure 12. NC discharge event. Upper line 
PMT, lower — NC 



unstability. Those periods may appear, for ex- 
ample, due to high humidity resulting in a high 
rate of breakdowns. Such periods are rejected if 
exist. Then a live time of measurement is calcu- 
lated taking into account a dead time, a rejected 
intervals and total number of frames recorded. At 
this stage a total number of normal single events 
in helium counters are fixed also in two energy re- 
gions — T nw in neutron window (35-180 ADC ch, 
or 170-870 keV) and T aw in a-window (181-248 
ADC ch, or 870-1200 keV). The T nw is corrected 
by a factor 1.1 of the neutron window accep- 
tance. The T aw is corrected to the interval 35-180 
ADC ch by a factor (248 - 181) / (180 - 35) = 2. 16. 
Both T nw and T aw are corrected by a factor 
19/(number of active helium counters). Also at 
this stage a relative uncertainty U es of those num- 
bers appearing due to uncertainty of the enegy 
scale is estimated. It is obtained as a variation of 
T nw and T aw when the boundaries of the windows 
were varied. 

Second, four PMT pulse height histograms are 
collected on the data being remained after rejec- 
tion. Initially, a threshold of a PMT signal to 
be accepted is obtained keeping in mind that the 
amount of double (triple, . . .) PMT events in a 



single frame should be less than (5-7)% of total 
PMT events. Usually it results in the threshold of 
10 ADC ch corresponding to 265 keV of electron 
scale, or 1.6 McV of neutron scale. Then, a si- 
multaneous events time window STW is obtained 
from the delay time distribution of PMT events 
coincident with NC events below the neutron win- 
dow (<170 keV). Usually STW = [-10; 6], or 
16 /its. Thus, one obtains an accidental ATW and 
a correlated CTW time windows also (usually 
CTW = [-160; -10] /us, ATW = [6; 160] fjs). Af- 
ter that, the PMT pulse height histograms PNC 
and PNA are collected for those PMT signals 
which are coincident with helium counter events 
in the neutron energy window and in the defi- 
nite time window (CTW and ATW respectively) . 
Similarly, the PAC and PAA are collected for 
those PMT signals which are coincident with 
helium counter events in the a-window and in 
the time windows CTW and ATW respectively. 
Again, the proper histograms are corrected by the 
neutron window acceptance, the number of ac- 
tive helium counters and the a-neutron windows 
width difference as well as it done for the T nw and 
T 

± aw • 

Third, the pure fast neutron signal is derived 
from the data. To make it clear one may look at 



10 









1 






100 150 



Figure 13. NC and PMT noise event. Upper line 
- PMT, lower — NC 



Figure 14. The model of distribution of time of 
delay in real measurement. See text for the labels 



the Figure El where a typical delay time distri- 
bution is simulated. It consists of pure fast neu- 
tron (region I) and pure Bi-Po (region II) events 
being coincident in the frame with specific time 
constants. Those events are grounded at a flat 
(with no characteristic time constant) distribu- 
tions of alphas (region III) and thermal neutrons 
(region IV) accidentally coincident with PMT sig- 
nals. So, the PNC histogram is a sum of pure fast 
neutrons, pure Bi-Po events and accidental back- 
ground (both a and thermal neutrons). The last 
component is obtained through direct measure- 
ment (the PNA histogram). The correlated Bi- 
Po histogram is obviously obtained as a remain- 
der of the PAC subtracted with the PAA. Thus, 
the fast neutron histogram FNS is obtained as 
follows: 



FNS = PNC - PNA - (PAC - PAA). 



(4) 



This procedure is applied to each particular en- 
ergy interval (see Table [2J of the spectrum. The 
uncertainties are statistical which propagate nor- 
mally and systematic one (U es , see above). The 
result of the procedure is corrected by the factor 
of the CTW window acceptance 1.38. Then the 
fast neutron flux is calculated with the efficiency 
averaged over each particular energy interval. 



Last, the slow neutron (i.e., under PMT thresh- 
old) signal is obtained as 



SNS = FNS - (Tnw - Taw), 



(5) 



each being corrected with proper efficiency. The 
FNS here is an integral of the PMT pulse height 
spectrum over all energy intervals. 

7. Example of data 

Let us consider a typical data measurement 
made at 660 m in low energy released region. 
Both PMT and NC channels were recorded in 
±100 mV dynamic range of the oscilloscope. This 
corresponds to [0.21-6.78] eeMeV (ee — electron 
equivalent) and [70-1240] keV ranges of the PMT 
and NC channels respectively. The frequency 
of digitizing was 3.125 kHz with 1024 points of 
frame (327.68 ^s long). The file J207020.osc was 
collected for 643 h during 2.07-29.07.2002 and 
contains 178285 frames. First analysis of the 
frame rate behavior revealed a period of unstabil- 
ity 336 — 297 = 39 h. Due to breakdowns at the 
period the frame rate increased drastically — see 
Figure ITK1 so the period was rejected. It resulted 
in 68085 frames remained after rejection. Each 
frame recording takes 60 ms, so the dead time 
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Figure 15. 3 He count rate time behavior at 660 m 



Figure 16. Pulse height spectrum of 3 He counters 
at 660 m 



appears to be 0.06 x 68085/3600 = 1.13 h. Thus, 
the live time is equal to 55 x (643-39) /60-1.13 = 
552.5 h. 

For NC normal single events the pulse height 
histogram is collected there fFigure lTBI) . One can 
clearly see the neutron signal expressed there as 
a typical peak at 760 keV. It is found that the 
histogram contains T nw = 17781 ± 133(stat) ± 
136(sysi) in the neutron window and T aw = 
6899 ± 83 {stat) ± 89 (syst) in the a- window. Af- 
ter corrections by a-to-neutron windows 2.164, 
by active neutron counters 19/18 = 1.055 and by 
neutron window acceptance 1.1 (only T nw ) and 
with quadratic sum of uncertainties one obtains 
T nw = 20646 ±221 and T aw = 15761 ±734 events. 

Next, the simultaneous time window STW 
has been checked for the correlated NC events 
in the energy region 70-170 keV (Figure 117(1 . 
One founds the STW = [-10; 6], or 16 fj,s, 
the CTW = [-160; -10] fis and the ATW = 
[6; 160] /is there. It's interesting to look at time 
delay distributions for the neutron window and 
a- window helium counters events (Figure^] and 
Figure 119(1 . Again, one can see clear fast neu- 
tron signature there. Then the PMT histograms 
were collected. Those histograms usually con- 
tain nothing specific — see, for example, Fig- 



ure 150] where the PMT pulse height spectrum of 
the events coincident with the NC pulses in the 
neutron window is presented. For particular re- 
gion (10-24 ADC ch, or 1.6-3.0 MeV of neutron 
energy) one finds after 1.055 correction applied 
the PNC = 572 ± 25 events as a correlated sig- 
nal. At the same region one finds the background 
PNA + (PAC - PA A) = 272 ± 32 events after 
the 1.055 and 2.164 corrections. Thus, the differ- 
ence corrected by 1.38 and 1.1 is a pure neutron 
signal FNS — 454 ± 62, which results in the neu- 
tron rate in this energy region of 0.82 ±0.11 h _1 . 
Then, using the efficiency 0.12 ±0.01 and the sur- 
face 6267 cm 2 one finds the flux in the energy re- 
gion of [3.04 ± 0.41] x 10~ 7 crn^s- 1 . The same 
procedure is applied to other energy regions. If a 
negative value is obtained there after subtraction 
one uses the uncertainty as an upper limit for the 
flux. 

Finally, one finds the slow neutron signal SNS. 
The integral of the fast neutron signal is FNS = 
734 ± 104 events. Thus, one finds the SNS = 
T nw - T aw - FNS = 4151 ± 319 events. With 
the efficiency 0.16 ± 0.03 it results in the flux 
[20.81 ± 0.60] x 10~ 7 cm- 2 s _1 in the neutron en- 
ergy region 0-1.6 MeV. The efficiency uncertainty 
isn't applied there. 
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Figure 17. The delay time distribution for the 
simultaneous window at 660 m 



8. Results and discussion 

An average rate of fast neutron detection events 
deep underground is something like ~1 h _1 giving 
thus several hundreds events per month roughly. 
Due to low statistics we perform result of mea- 
surements in different energy ranges rather than 
spectrum plots. 

Table 0] contains results of the measurement 
of the neutron background at different levels of 
the Pyhasalmi mine. Uncertainties expressed ex- 
plicitely includes both statistical and systemati- 
cal ones. In addition, there are uncertainties of 
average detection efficiency for every energy re- 
gion. The neutron energy above ~1.5 MeV has 
been measured directly as an amplitude of light 
flash in the scintillator coincident with capture in 
3 He counters. Count rate of the coincident events 
is a base for flux calculation. Neutron with en- 
ergy <1.5 MeV does not produce detectable light 
flash but can be effectively detected by helium 
counters. The neutron flux with those energies is 
deduced from count rate of helium counters sub- 
tracted by fast neutron rate. Those fluxes that 
appeared to be under uncertainties are expressed 
as an upper limits. 

At 400 m level one can see significant excess 
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Figure 18. The delay time distribution for the 
neutron window at 660 m 



in the energy region >25 MeV. An obvious rea- 
son is a production of neutrons in hadron showers 
generated by cosmic muons in the lead shielding. 
Muon induced neutrons mostly have an evaporat- 
ing spectrum that is their energy lies in the range 
1-10 MeV, and in principle can be measured in 
the detector. But in this case the energy released 
by evaporated neutrons in the scintillator is ac- 
companied with and usually much less than the 
energy deposited by muon induced hadron show- 
ers. So, it can't be distinguished. Nevertheless, 
we put those neutrons in high energy released re- 
gion just to separate from other ones. 

Significant excess of neutrons at 1410 m can 
be easily reasoned by higher contamination of 
U/Th in the decorative granite powder covering 
the walls of the room where the measurement was 
done. 
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Table 4 



Neutron background at different levels of the Pyhasalmi mine, 10 7 cm 2 s 
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Energy released region 
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0.12±0.01 


0.10±0.01 
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Figure 19. The delay time distribution for the 
a-window at 660 m 



Figure 20. The pulse height histogram of the 
PMT events correlated with the NC events in the 
neutron window at 660 m 
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